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OCHPARISOH OF CALCULATSD AMD KSASURED 
KCDSL ROTOR LOADING AHD WAKE GEOMETRY 


Wayne Johnson* 

Aaeo Research Center 
and 

Aeromechanics Laboratory 
AVRADCOM Research and Technology Laboratories 

ABSTRACT 

Calculated blade bound circulation and woke geomotry arc compared 
,ith measured results for a model helicopter rotor in hover and forward 
flight. Hover results are presented for rectangular tip and ogee tip planform 
blades. The correlation is quite good uhen the meacured wake geometry 
chaiacteristics are used in the analysts. Available prescribed uakc 
geometry models are found to give fair predictions of the loading, but 
they do not produce a reasonable prediction of tho induced pouer. Fon^arcl 
•flight results arc presented for twisted and untwisted blades. Fair 
correlation between measurements and calculations is found for the bound 
circulation distribution on the advancing side. Tl»e tip vortex geometry 
in the vicinity of the advancing blade in forward flight was predicted 
well by the free wake calculation used, although th-e wake geomotry did not 
have a significant influence on the calculated loading and psrforranoe for 
"tho cD.r3GS CO ns id-Grcrt • 


*Head of Rotorcraft Rer.earch 


Section, Largc-Scalo Aerodynamics Branch 



Ih’TOOKJCTION 


A series of oxperlaento have been conducted at Aess Research 

Center, In which a laser voloclnotor was eaployed to ncaouro tha flow 

field of a Rodol helicopter rotor. Tno results of these exporinents have 

included neasured bound circulation and tip vortex location for a hoverins 

12 

rotor with two blade tip planfoms; * noaoured bound circulation on the advancln 

side of the rotor in forward fllghtr and measured tip vertex position on 

4 5 

the advancing side in fortiard flight, for two blado twist toIuco. 

Each of these experiments used a two-bladed rotor with a radius of approximately 
1 meter, operated at low tip speed. The forward flight exporimonts were 
conducted in a 7- by 10-ft wind tunnel. A two-coraponent laser vcloclmotcr 
was used to ncacure the flovr field in the rotor wake. Tho blade loading 
was obtained by measuring the circulation around a box enclosing tho blade 
at a spocified radial station. Tho wake geometry, in toms of the location 
of the tip vortices, was obtained from measurements of the induced velocity 
along radial and vertical traverses tteough the wake. Ccnplctc details of 
the experiments are given in the original papers.^ ^ 

A principal objective of this scries of experinonts has been to 
provide data for the development and verification of rotor analytical 
models. The purpose of the present paper is to compare the measured 
loading and wake geometry vrith calculated results. For tho calculations, 
a recently developed helicopter analysts was used.^ The rotor wake analysis 
was based on a discrete element representation of the vorticity, vrith modolo 
for the wake rollup and the distorted wake geometry. The loading was calculated 
using lifting lino theory, vrith corrections available for three-dimensional 
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effects at the tip and at vortex/t>lsdo Intcractiono. The rotor blade 

flappl»JU notion was calculated by means of a harmonic analyclo nethod. 

An exlctins analyclo was used to calculate tho distorted, free vfsko 

7 

geoBotry In forward flight.' 

AHALYTIC/\L tiOISL 

A top view of the waho nodel in shown In figure 1. Tho Induced 
velocity, circulation, and lift aro evaluated at 15 stations along the 
blade, as shown by the dots. The trailed vortlclty in the ifake directly 
behind, the blado is represented by discrete vortex lines positioned nldway 
between tho points at which the circulation is calculated. The strength 
of each trailed lineals defined by the difference between the bound 
circulation at successive radial stations. This jart of the Eodel Is a 
coBiJion numerical implementation of lifting line theory. For tho rotary 
wing it is necessary also to model tVve rollup of the vortlclty into a 
concentrated tip vortex, because the blade encounters tho vrakc from tho 
preceding blades as it rotates. Therefore, after an asimuthal extent of 
30° in this case, tho trailed vortlclty is concentrated into a slriglc tip 
vortex lino,- with strength equal to the Raximim bound circulation of the 
blade. The azimuthal extent of the near vako was varied from 15° to 60° 

(and the number of radial stations from 10 to 15) with little effect on tho 
calculated loading. Tho key element In the model is the vortlclty distribution 
in the wake when it encounters the following blade, at which point it must 
be concentrated into the tip vortex. For conservation of vorticity, thoro 
must bo an Inboarrl sheet of trailed vorticity vilth equal total strength aM 
opposite sign as the tip vortex. This inboard sheet is much less important 
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than tha tip vortex, partly because tho vortlclty la not no concentrated 

nrd jartly because it la convoctod doirnward at about twice the rate as the 
8 

tip vortex, Conocqucntly the entire inboard sheet can bo reasonably represfinted 
ao a single vortex lino, shown dashed in figuro 1, with a largo core radius 
of the rotor radius in this case) to avoid unrealistically large 
induced velocities when the waits passes under tho following blade, A 
similar wake model is used for the other blade of the rotor, Tho detailed 
near wako model need not bo used for the second blade since it is far from 
whore the induced velocity is being calculated on the first blade (figure l). 

Tho curved trailed vortex elements in the waltc are represented by a connected 
series of straight line segments, with an azimuthal increment of 15° in this 
case. Five revolutions of tho walte behind each blade are modelled In this 
fashion (only the first one-half revolution is shown in figure l). An 
additional 30 revolutions of the t?ako are modelled using rectangular vortex 
shoet panels to construct a cylinder of axial and circumferential vorticity 
representing tha tip vortices, and an axial line vortex representing 
the root vortices. By this moans the wsjee far from the rotor blade (extending 
approximately from 2, to 10 rotor i-adii below the disk) is economically 
accounted for; if this part of the wake were neglected, the induced velocity 
at the rotor disk would bo significantly underestimated, 

Tho geomeiry of tho tip vortex in hover is described by the following 
coordinates » 

X = (it D^) cos('^i>-4) 

y = (1 + D^) sln(‘f>- (^) 

z « 
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A tip-path piano cooirdlnato system vlth origin centered on tho rotor shaft 
is used I X and y aro rectangular coordinates in tho rotor piano t s io 


measured vertically upward from tho piano of the blado tijxi* The aairiuth 

angle of tho rotor blado io , measured from tho x axlej ^ is the 

Make age, i.o. the azimuth angle along the valto holix, measured from the 

rotor blade# D and D are the vertical and radial distortions of tho uake 
2 r 

duo to its self- induced velocity. A two-stage vortical convection and 
exponential radial contraction model is used for hovori 



r -Kjf 

\ -KjTr 

-(1 - e-%‘^)(l - K^) 




(a similar model is used for the distortion of tho inside and outside 
edges of ths Inboard wake sheet, but only tho tip vortox geometry has a 
significant effect on the loading). Ths parameters K^, K^, and dotemine 
tho position of the tip vortex whan it encounters the follouing blade? 
determines the contracted radius in ths far wake; and K„ determines tho 

t 

vertical convection after the first blade pjassage. An undlstoi?,cd wake geometry 
can be obUinod using *= 0 and D^, = - V . vfhere Is the mean 

inflow ratio at the rotor disk. 


A similar walce model is used for the forward flight calculations. 
The bound circulation in forward flight varies azlmuthally as well as 
radially, so the shod wake Is also included. Tho shed wake is modelled by 
radial, dlsoreto vortox lines. The vortlclty ctrensth varies linearly 
along each vortex line element, both trailed and shed. In forviard flight 
it is only necessary to use two revolutions of ths wake behind each blade 
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to ealciOato ths induced velocity. Tho scene try of tho tip vortex rolGtive 
to tho tip-rath piano lo described by tl» followlns coordL-Btcoi 
«tpp“ *Dy(K\f) 

^tpp “ + Dy(^^» .‘j5>) 

=tpp“>/^*^n%p +Dg(fJ^«S>) 

Hhoro ^ is tho rotor advance ratio (for,;ard cr< 5 od divided by tip epood) 

**'tpp of attack of tho tlp-path piano, positlvo for 

rcarifard tilt. Tho coordinate s is positive dounstrean, and y* is 
positive to tho richt. Hind or tunnsl axis coordinates aro obtained by 
accountins for the blade flapping notions 

»hore (Ej, Is the bleat, itrocono assle, e„s lonsltadlsel 

and lateral flapping angles (relative to the shaft), and Is tho shaft 
angle of attack. An undistorted yolsc gconotry, that only accounts for the 
vertical convection due to the ncan induced velocity, is obtained 



0 and D, 



A free wake geonetxy calculation can also 


ba uGGd. 


■ ■ .jjovEa RESULTS 

Exporimcntal data are available for tho circulation and wake geometry 
In hover with two tip sha-pos, rectangular and ogee.^’^ Tho principal 
parameters defining tho rotor and oporatlng conditions aro given in Table if 
further details, particularly for tho twist dlstrlbuUoa and ogoo tip planform 
are given in the original papers. Tho ogee tip planforn extends over the 
outer IQ?-; span of the blade. The solidity used for the ogee tip is tae 
thrust-equivalent values (radius -squared weighted). In all calculations, 
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Table 1. Rotor paracotors and operating condition for hover caoes. 



Rectangular 
Tip Blades 

Ogee Tip 
Blades 

Blade radius, R 

1,045 la 

1.045 la 

Rotor tip speed 

76.6 n/sec 

76.6 e/bcc 

Number of blades 

2 

2 

Chord, c/R 

0.0729 

0.0729 

Precono 

1 . 5 ° 

1 . 5 ° 

Twist 

0 

1 

0 

1 

Solidity, 

0.0464 

0.0396 

Thrust, 

0.099 

0.103 
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tho rotor thrist was trlnaod to ths value given In 'teble It which was 
obtained by integrating tho measured loading distribution. 

Tte calculated and measured circulation distributions are compared 
in figure 2 for the rectangular tip blade in hover. The calculations 
used tho measured position of the tip vortex vihen it first passed under tho 
following blado, and a value of Kg “ 0.0550. A tip loss factor of B « 0.985 
was used. The correlation exhibited in figure 2 is quite good. The wake 
geometry used is shown in figure 3. in terras of the intersections of the 
tip vortices with a vertical plane under the rotor blado (the x-y plans, 
with tho blade at azimuth angle 90°). The undlstorted iralio geometry 
is also shown in figure 3, and figure 4 shows the corresponding calculated 
circulation. Without a propor modelling of the walte contraction, the 
calculated loading distribution is very different from the measured distribution 

Prescribed vraltc geometry models for hovering rotors havo been 
developed by Landgrobe^ and by Kocurek and Tangier.'® They developed 
enpirlcal oxprossions for the constants K^, K^, and by matching 
rotor wal:e geometry data obtained using flow visualization techniques. 

Figures 5 and 6 show the loading calculated using these two wake goonetry 
models. For the present case, the only significant difference between 
the Landgrebs model and the Kocurek and Tangier model is in the jaramotcr K^, 
hence in the vertical separation of the vortex when it first encounters the 
following blado (see figure 3)* Th° Landgrebs model places tho vortex a 
small distance R) closer to tho blade, hence the peak bound circulation 

is somovjhat larger tton for the Kocurek and Tangier model (compare figuie^ 
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5 and 6)* Thsso prescri'bsd tfaI;o nodelo give a valua of Kg •• 0t0608» 

Figure 7 shows the loading calculated uoing tho iscasurcd location for tha 

tip vortcK when it first piseeQ under the blade, and a value of Kg " 0#0608« 

Figures 5, 6, and ? differ only in the location of tho first vortex* under 

9 10 

the blade. The data used to construct the prescribed wake models, ’ - 

show a scatter In tho vortex position of perhaps si 1 to 2^ of the blade 
radius. Within that accuracy, tho tip vortex geonetrlos defined by the 
two models aro Identical for this case, and both correctly predict the 
location of tho first vortex under the blade. However, although tho difference! 
in the vortex positions are slight, using tho measured position (figure ?) 
definitely Improved tho calculated loading distribution, compared with using 
the prescribed wake models (figures 5 and 6). 

Figures 2 and 7 differ only in ths value of tho parameter Kg, 

Which Is the vertical convection rate aftor the vortex first encounters 
the following blade . Tno measured urWc geonotry data does not provide a 
value for Kg. Initially therefore the loading t:as calculated using tho 
same value of as given by tho prescribed viakc models (figure ?). By 
using a iO^ smaller value of K^, . tho correlation of the measured and 
calculated circulation distributions was defitiitely Iminroved (figure 2). 

The calculated performance provides further support for tho use of a smaller 
value of Ko than obtained from the prescribod tfake modelG. Using K„ « O.C608 
(figure ?), tho calculated Induced pov'or was only 3 % higher than the ideal 
momentum- theory Induced loss, which is much bettor parformanco than would 
be expected, especially with two blades. Using K 2 “ O. 055 O (figure 2), tho 
calculated induced power vfas higher than Ideal, a more reasonablo value, 
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The calculated and nbzusurod circulation distributions arc comiared 
in figure 8 for the ogee tip blade in hover. The calculations used the 
rieasurcd position of tho tip vortex when it first passed under the following 
blade, and a value of Kg •» O.O5OO (figure 3). With the ogee planforn, the 
measured data show the tip vortex rollup occurring at a radial station of 
R, It was possible to proixsrly model this rollup position in the 
analysis by simply using a tip loss factor of B “ 0.9^. . Tto correlation 
exhibited in figure 8 is nearly as good as for tho rectangular tip blade. 

For the ogee tip there aare no corresponding prescribed wake models to provide 
a guide for the choice of K^; the wake convection must oa specified on the 
basis of tho peak circulation and the predicted induced power. Properly 
modelling tho tip vortex rollup at R is essential for an accurate 
calculation of the loading distribution. For comparison, figure 9 shows the 
loading calculated for the ogee tip when the vortex rollup renains at r/R “ 1 
as for th.e rectangular blade. It is interesting to note that the ogee tip 
blade can be considered equivalent to a rectangular tip blade with a radius - 
of R as regards the bound circulation distribution and uaks geometry;^ 
and as regards tiic rotor performance parameters (disk loading, induced power, 
and figure of merit). Of course, the egos planform in addition produces 
a tip vortex with a larger core radius and smaller peak velocity. 

FORWARD FLIGHT RESULTS 

Experimental data are available for tho circulation and wake 
geometry in forward flight for blades with ssro twist and -11° twist. 

The principal parameters defining the rotor and operating conditions are 
given in Table 2. Tho measured bound circulation is available only for tho 
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TStle 2. Rotor paronotoro ood oporatlng condition tor toronrd tllcM cases. 


Collective pitch, 

Tip- path plane angle, 



Zero Twist 

-11 THISX 


Blades 

Blades 

Blade radius, B 

1 .066 m 

1.045 n 

Rotor tip speed 

6?.0 m/sec 

65.7 n/sec 

Nuiaber of blades 

2 

• 2 

Chord, c/R 

0.1012 

0.0729 

Precone 

1.5° 

1.5° 

Solidity,^ 

0.0644 

0.0464 

Lock nviraber, % 

3.3 

3.78 

Advance ratio, V/X^R 

0.18 

0.18 

Shaft angle, 

-id° 

-10° 

r\ 


8.5 

- 6 . 6 ° 


10. r 

-b-f 
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untsflGtod \>lado, on tha advanolns oido ot 90°. Tho ncasurcd end calculated 
circulation distributions exo coapaxed in figur© 10 • Ths circulation box 
encloses sons of tho shod traI;o vorticity, which is otiui valent to measuring 
tha circulation of the blado iihsn it ia at a alightly earlier osieuth position; 
hence the calculated loading at “ 75° is also Bhotm. Tho correlation 
exhibited in figure 10 is fair, although both tho ccatter in tho neaGured 
data and tha differences between tho moaoured and calculated results are 
larger than for ths hovor eases, Ths calculated loading in figure 10 was 
obtained using the free wake geometry model, but the results obtained 
using an undistortod wake geometry were nearly tha sane. Using tho two 
wake models, only small differences were produced in the calculated 
tlurust, power, ard flapping, Kowover, using nonunlfom wake- induced 
velocity produced largo dii'feroncss compared with using uniform tnflovj 
(with uniform induced velocity, tho bound circulation at = 90 ® would bo 
a linear function of tho radial otation r/n). 

Keasurements arc available for the geometry of the tip vortex from 
tho blade at 270®, in tha vicinity of the other blade at V 90® 

on the advancing slde.^ Figures 11 and 12 compare the calculated free va!:o 
geometry results with the measured tip vortex position for tho untwisted 
blades. A top view is shown in figure 11 and a side view in figure 12. 

A tunnel or wind axis system is used, with the x coordinate directed 
downstream, the y coor^Unate toward the advancing side, and tho s coordinate 
upward. The rotor hub Is the center of tho coordinate cyctem. Tlw waI<o 
geometry is calculated at 15® azli.iuth intort'als, so the range of the experimental 
data is covered by just two vortex lino aegments. The corresponding results 
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for tho -11“ t»lct M»dco oro stora In flguroo 13 nnd 14. Flsnio. 12 an,l 14 
alco nlB« tte ecoBCtr, obtained unlns ««. undlotortod «nl=0 note!. Thn 
correlation botnoon tho nonourod and ealculalod free anUo eooaetry lo Good. ^ 
conaldorlns that the calcalaUon doca not attonpt to nodol th) flno dotollc 
or the eoonotry. At thla point la the Kako. tho predicted dletortlon la 
onaU in tho « and y dlroctlonn. Tho vorlleal convoetlon lo nlenltleant 
hooovor. and has tho not offoct oC kooplns tho vortex In tho vlelnlty of tho 
blade at 'f- 90° nearly parallel to tho tlp-path piano (nhlch la tilted 
fervard rolatlvo to tho tunnel area). Thin trend Is conflmod by the 
aeasured pooltlono. In eontiuat. the undlotortod Goonotry shoun tho tip 
vortex bcln« steadily eoaveetod dounvard rolatlvo to the tlp-pnth piano as 


tho element age Increanen. 

It oust bo recoonlood houovcr that tlie exporlnontol data doos not 
provide a partloularly severe tent of tho free u=P.o caleulatlon. The tip-path 
plane ansle of attae!. Is fairly lardO , anl tho tip vortex oneouatera tho 
followla-, blr.de ultl. nub.-,tantl.al vertical separation. Conoe,u=ntly tho tashs 
of caloulatluG tho uako G«o=olxy and lo,adlnG are easier than for snail tlp-path 
plane anGle of atUek. In fact. althoUEh tho use of nouualforn Inflou lo 
inportant. there Is little dlfforonoe for this ease botueen the free v.ako 
results and the undlstorted uaho results, as rosards the porforaanee ,an.l 

blndc loading* 


COllCLUSICNS AND RECOW'IEHD.\TIONS 
Ce-ilculatcd blade bound circulation and vako geonetry have been 


con pared ulth neac>uret^ 


renult-n for a model holicortor rotor in hover and 
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for»;ard fllGht. Quito good corrolation wao found for the predlctud and 
noasurcd circulation diotrl'butioiKJ in hovor. Thlo corrolation was only 
pocciblo using tha neaaurcd waI;o gcoaotry charactcrlntlcs, including the tip 
vorten rollup, the radial and vortical position of tho vortex when it 
first encounters tho following blade, and tha vortical convoctlon rate in 
the far wako. An accurato vralte goonetry calculation will ba necessary for a 
coiapletely analytical calculation of hover perfonnanco and loading. The 
waJ;e goonetry defined by tho available prescribed wake nodels was found to 
give tha general features of tho loading distribution well, although significant 
cianges in tha paak loading are produced by georaetry differences that are 
within tho scatter of tha data from which the prescribed wake nod ols were 
constructed. Moreover, a reasonable prediction of the induced power (which 
is veri* sensitive to the ifako goonetry changes) was not obtained tislng the 
prescribed walcc nodels. The correlation between tho calculated and 
neasured circulation distributions In forward flight was acceptable, but 
not as good as for Uie hover cases. Tlie tip vortex gconetry in the vicinity 
of the advancing blade in foniard flight vras predicted fairly well by the 
free wake calculation used,' although for these cases the wake geonetry diid 
not have a significant influence on the calculated loading and perfomanoe. 

The experimental program that produced the data examined in this paper 
Is continuing. The work in hover will be expanded to include variations in 
the number of blades, thrust, planform, and tip Mach namber. It should be 
pooolblo in these tests to also measure the vertical convection rate after the 
tip vortex first encounters the following blade. In future forward flight 
tests, it villi be possible to measure the loading at azimuthal stations on 
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the front and roar of tho dicki by uains a three-component laser velocioetcr 
systen. It is rccoRncndcd that data bs obtalnod at only a few azimuth ctationa, 
for a wide reside of rotor oporating coiditlons (as opposed to collecting 
deta at enough azimuth stations to define the loading distribution over the 
entire disk, but only for one or tr;o operating conditions)* Even more useful 
would be the loading at a fixed radial station as a function of azimuth, since 
vortex- induced loads are noiro readily discerned when the loading Is viewed 
as a function of azimuth rather than radius. The operating conditions 
should Include cases of small tip- pith plane incidence, to provide a better 
test of the loading and wake geometry predictions. 
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Figuro 2« Conparison of neasurcd and calculated circulation for rectangular 


tip blade in hover, Kith meacurcd waJce geometry and « 0.0550 
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Figure 3. Location of tip vortices in hover (intersections with vertical 

plane under rotor blade? flagged points obtained by LV neasurenents) 
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au/J 



Figure 7. 


Bound circulation calculated 


vjith neaeured 


wahe 


geometry 



Figxire 8 



Conparison of neasured and calculated circulation for ogee tip 


blade in hover 




r/R 


Figure 9. Bound circulation calculated for ogee tip Blade wltVi tip vortex 
rollup at r/R =• 1. 



au/J 



Figure 10. Gonparir.on of tnear.ured and calculated circulation on the advancing 
side in forward flight (v/aR = O.lR, untwisted blade) 
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O MEASURED 

FREE WAKE CALCULATION 

— — UNDISTORTED WAKE 



FiGure 12. Side view of tip vertex geometry in vicinity of blade at 90° 

(v/iTLR = 0.18, untwisted blade) 
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O TREASURED 

FREE WAKE CALCULATION 

UNDISTORTED WAKE 


BLADE 

CHORD 



Figure l4. Side view of tip vortex geometry in vicinity of blade at ^'■p = 90° 
(V/n.R = 0ll8, -11° ’ twist blade) ; 


- 31 - 







.l^NGLEY RESE/ 



76 00 


88 9444 



